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Abstract. We develop a general theory of "almost Hadamard matrices". These are by 
definition the matrices H S Mn{R) having the property that U = H/^/N is orthogonal, 
and is a local maximum of the 1-norm on 0{N). Our study includes a detailed discussion 
of the circulant case (iJ^ = Ij-i) of the two-entry case {Hij G {x,y}), with the 
construction of several families of examples, and some 1-norm computations. 



Introduction 

An Hadamard matrix is a square matrix having ±1 entries, whose rows are pairwise 
orthogonal The simplest example, appearing at = 2, is the Walsh matrix: 



Ho 



1 1 
1 -1 



At = 3 we cannot have examples, due to the orthogonality condition, which forces 
to be even. At = 4 now, we have several examples, for instance = H2® H2'- 



Ha 



(\ 
1 
1 



1 \ 

-1 
-1 

1/ 



For higher values of A^, the construction of Hadamard matrices is quite a tricky problem. 
First, by permuting rows and columns or by multiplying them by —1, we can always 
assume that the first 3 rows of our matrix look as follows: 

1 

1 -1 . 

- 1 1 



H 



1. 
1. 



.1 
.1 



-1 



,1 1 
- 1 -1 . 
,1 -1. 



.1 \ 

- 1 

- 1 



v 



Now if we denote by x, y, z, t the sizes of the 4 columns, the orthogonality conditions 
between the first 3 rows give x = y = z = t^'s,oN = x^y^z^tv&& multiple of 4. 
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A similar analysis with 4 rows instead of 3, or any other kind of abstract or concrete 
consideration doesn't give any further restriction on N, and we have: 

Hadamard Conjecture (HC). There is at least one Hadamard matrix of size N x N, 
for any N e AN. 

This conjecture, going back to the 19th century, is probably one of the most beautiful 
statements in combinatorics, and in mathematics in general. The numeric verification so 
far goes up to = 664, see [24]. For a general presentation of the problem, see ^5] . 

At the level of concrete examples, the only ones which are simple to describe are the 
tensor powers of the Walsh matrix, having size 2^^. For some other examples, see [20] . 

Yet another good problem, simple-looking as well, concerns the circulant case. Given 
a vector 7 G (±1)^, one can ask whether the matrix H G M7v(±l) defined by Hij = 7j_j 
is Hadamard or not. Here is a solution to the problem, appearing at A^ = 4: 



/-I 1 


1 1\ 


1 -1 


1 1 


1 1 


-1 1 


u 1 


1 -V 



More generally, any vector 7 G (±1)^ satisfying Ylli — ^ solution to the problem. 
The following conjecture, going back to [M], states that there are no other solutions: 

Circulant Hadamard Conjecture (CHC). There is no circulant Hadamard matrix of 
size N X N, for any N A. 

The fact that such a simple-looking problem is still open might seem quite surprizing. 
If we denote byS'c{l,...,A^} the set of positions of the —1 entries of 7, the Hadamard 
matrix condition is simply jS" fl (S* -|- fc)| = 15*1 — A^/4, for any k 0, taken modulo A^. 
Thus, the above conjecture simply states that at A^ 7^ 4, such a set S cannot exist! 

Summarizing, the Hadamard matrices are very easy to introduce, and they quickly 
lead to some difficult and interesting combinatorial problems. Regarding now the other 
motivations for studying such matrices, these are quite varied: 

(1) The Hadamard matrices were first studied by Sylvester [31], who was seemingly 
attracted by their plastic beauty: just replace the ±1 entries by black and white 
tiles, and admire the symmetries and dissymmetries of the resulting design! 

(2) More concretely now, the Hadamard matrices can be used for various coding pur- 
poses, and have several applications to engineering, and quantum physics. For 
instance the Walsh matrices H2'' are used in the Olivia MFSK radio protocol. 

Most applications of the Hadamard matrices, however, come from their generalizations. 
A "complex Hadamard" matrix is a matrix H G Mn{C) all whose entries are on the unit 
circle, and whose rows are pairwise orthogonal. The basic example is F/v = y/NFj^, where 
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Fn is the matrix of the Fourier transform over Zat. That is, with cu = e^'^*/^: 



/I 


1 


1 


•• 1 \ 


1 


CO 






1 

















As a first observation, the existence of this matrix prevents the existence of a "com- 
plex version" of the HC. However, when trying to construct complex Hadamard matrices 
by using roots of unity of a given order, a wide, subtle, and quite poorly understood 
generalization of the HC problematics appears. See e.g. [12], [26], [27] . 

As for the motivations and applications, these partly come from pure mathematics, cf. 
e.g. [S], [IZ], [23], [32], and partly come from quantum physics, cf. e.g. [S], [3Z]- 

Let us go back now to the real case. Since the determinant of vectors is maximized 
when these vectors are chosen pairwise orthogonal, we have the following result: 

Theorem A. For a matrix H E MAr(±l) we have \ detif| < N^^"^, with equality if and 
only if H is Hadamard. 

This result, due to Hadamard himself [12], has led to a number of interesting problems, 
and to the general development of the theory of Hadamard matrices. See |2U] . 

As already mentioned, in order for an Hadamard matrix to exist, its size A^ must 
be a multiple of 4. For numbers of type N = An + k with k = 1,2, 3, several "real" 
generalizations of the Hadamard matrices have been constructed. The idea is usually to 
consider matrices H e M7v(±l), whose rows are as orthogonal as they can be: 

Definition A. A "quasi- Hadamard" matrix is a square matrix H G M7v(±l) which is as 
orthogonal as possible, e.g. which maximizes the quantity \detH\. 

This definition is of course a bit vague, but the main idea is there. For a detailed 
discussion of the different notions here, we refer to the articles [1], [25], [31] . 

Yet another interpretation of the Hadamard matrices comes from the Cauchy-Schwarz 
inequality. Since for U G 0{N) we have \ \U\\2 = VN, we obtain: 

Theorem B. For a matrix U G 0{N) we have \ \U\\i < N\/N, with equality if and only 
if H = y/NU is Hadamard. 

This result, first pointed out in [5], shows that the matrices of type H = y/NU , with 
U G 0{N) being a maximizer of the 1-norm on 0{N), can be thought of as being some 
kind of "analytic generalizations" of the Hadamard matrices. Note that such matrices 
exist for any A^, in particular for N = 4n + k with k = 1,2, 3. 
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As an example, the maximum of the 1-norm on 0(3) can be shown to be the number 
5, coming from the matrix Us = K^/y/S and its various conjugates, where: 




This resuh, proved in |5] by using the Euler-Rodrigues formula, is actually something 
quite accidental. In general, the integration on 0{N) is quite a subtle business, and the 
maximum of the 1-norm is quite difficult to approach. See [6], [7], |15]. 

These integration problems make the above type of matrices quite hard to investigate. 
Instead of looking directly at them, we will rather enlarge the attention to the matrices of 
type H = \/NU, where U G 0{N) is a local maximizer of the 1-norm on 0{N). Indeed, 
according to the Hessian computations in [S], these latter matrices are characterized by 



the fact that all their entries are nonzero, and SU^ > 0, where Sij = sgn(t/, 



Summarizing, Theorem B suggests the following definition: 



Definition B. An "almost Hadamard" matrix is a square matrix H G Mjv(M) such that 
U = H/^/N is orthogonal, and is a local maximum of the 1-norm on 0{N). 



The basic examples are of course the Hadamard matrices. We have as well the following 
N X N matrix, with iV G N arbitrary, which generalizes the above matrices K^, K^. 



(2-N 2 
2 2 - 



K 



N 



\ 



\ 



N) 



A lot of other interesting examples exist, as we will show in this paper. Here is for 
instance a remarkable matrix, having order A^ G 2N + 1, and circulant structure: 



COS" 



^1 {N-l)-K 
N 



COS 



-1 Tl_ 

N 



— COS 

1 

COS" 



-1 Tl_ 

N 



-1 27r 

N 



COS 

— COS 

— COS" 



1 27r 

N 

-1 TT 



N 



-1 Stt 

N 



COS 
— COS 



-1 (jV-l)7r \ 
N 

1 {N-2)-K 



N 



Yet another series, with A^ = -|- g -|- 1, where q = is a. prime power, comes from 
the adjacency matrix of the projective plane over F^. Here is for instance the matrix 



ALMOST HADAMARD MATRICES 



5 



associated to the Fano plane {q — 2), where x = 2 — 4\/2, y = 2 + 3\/2: 

y y y y\ 

X X y y 
y y X X 
X y y X 
y X X y 
y X y X 
X y X yj 

The aim of the present paper is to provide a systematic study of such matrices, with 
the construction of a number of non-trivial examples, and with the development of some 
general theory as well. Our motivation comes from two kinds of problems: 

(1) The Hadamard matrix problematics. The world of Hadamard matrices is ex- 
tremely rigid, and we think that our study of almost Hadamard matrices, where 
there is much more freedom, can help. As an example, there are several non- 
trivial classes of circulant almost Hadamard matrices at any G N. and trying 
to understand them might end up in sheding some new light on the CHC. 

(2) Generalizations of Hadamard matrices. The Hadamard matrices have applications 
in a number of areas of physics and engineering, notably in coding theory, and in 
various branches of quantum physics. One problem, however, is that these matrices 
exist only at iV = An. At N = An + k with A; = 1, 2, 3 some generalizations would 
be needed, and we believe that our almost Hadamard matrices can help. 

The paper is organized as follows: 1 is a preliminary section, in 2-3 we investigate two 
special cases, namely the circulant case and the two-entry case, and 4 contains a list of 
examples. The final section, 5, contains a few concluding remarks. 

Acknowledgements. We would like to thank Guillaume Aubrun for several useful dis- 
cussions. The work of T.B. was supported by the ANR grant "Granma". I.N. acknowl- 
edges financial support from the ANR project OSvsQPI 2011 BSOl 008 01 and from 
a CNRS PEPS grant. The work of K.Z. was supported by the grant N N202 261938, 
financed by the Polish Ministry of Science and Higher Education. 

1. Preliminaries 

We consider in this paper various square matrices M e Miv(C). The indices of our 
matrices will usually range in the set {0, 1, . . . , A" — 1}. 

Definition 1.1. We use the following special N x N matrices: 

(1) Iat is the identity matrix. 

(2) Jjv is the "flat" matrix, having all entries equal to 1/N. 

(3) Fjv is the Fourier matrix, given by (i^Ar)ij = u'^^ /y/N, with u — e^'^'^l^ . 



h 



2^/^ 



fx 


X 


X 


X 


y 


y 


X 


y 


y 


y 


X 


y 


y 


X 


y 


y 


y 


X 


\y 


y 


X 
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We denote by D the generic diagonal matrices, by U the generic orthogonal or unitary 
matrices, and by H the generic Hadamard matrices, and their generalizations. 
Our starting point is the following observation, from [5]: 

Proposition 1.2. For U G 0{N) we have \\U\\i < N\fN , with equality if and only if 
H = \/NU is Hadamard. 

Proof. The first assertion follows from the Cauchy-Schwarz inequality: 

ij \ ij / \ ij / 

For having equality the numbers \Uij\ must be equal, and since the sum of squares of 
these numbers is A^, we must have \Uij\ = 1/\/N, which gives the result. □ 

As already mentioned in the introduction, the study of maximizers of the 1-norm on 
0{N) is a quite difficult task. However, we have here the following result, from [5]: 

Theorem 1.3. For U G 0{N), the following are equivalent: 

(1) U is a local maximizer of the 1-norm on 0{N). 

(2) Uij 7^ 0, and SU* > 0, where Sij = sgn(f/jj). 

Proof. As already mentioned, this result is from [5]. Here is the idea of the proof: 

Let us first prove that if ?7 is a local maximizer of the 1-norm, then Uij ^ 0. For this 
purpose, let Ui, . . . , Un be the columns of U, and let us perform a rotation of f/i, f/2: 

f Ul\ _ /cost ■ f/i - sint • 
\U^J ~ \smt- Ui + cost ■ U2) 

In order to compute the 1-norm, let us permute the columns of U, in such a way that 
the first two rows look as follows, with Xk 7^ 0, Yk 7^ 0, AkCk > 0, B^Dk < 0: 

fUi\ _ fO Y A B\ 
\U2) ~ [0 X C DJ 

If we agree that the lower-case letters denote the 1-norms of the corresponding upper- 
case vectors, and we let K = U3 + ... + un , then for t > small we have: 

||f/*||i = II cost ■ f/i - sint ■ f/2||i + II sint ■ f/i + cost ■ f/2||i + 

= (cos t + smt){x + y + a + d) + (cos t — sin t){b + c) + K 

Now since U locally maximizes the 1-norm on 0{N), the derivative of this quantity 
must be negative in the limit t — t- 0. So, we obtain the following inequality: 

{x + y + a + d)-{b + c)<0 
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Consider now the matrix obtained by interchanging Ui, U2- Since this matrix must be 
as well a local maximizer of the 1-norm on 0{N), we obtain: 

(x + y + b + c)-{a + d) <0 

We deduce that x + y — 0, so x — y = 0, and the entries of Ui, U2 must appear at the 
same positions. By permuting the rows of U the same must hold for any two rows Ui, Uj. 

Now since U cannot have zero columns, all its entries must be nonzero, as claimed. 

It remains to prove that if Uij 7^ 0, then U is a local maximizer of F{U) = \ \U\\i if and 
only if SU^ > 0, where Sij = sgn{Uij). For this purpose, we differentiate F: 

dF — SijdUij 

We know that 0{N) consists of the zeroes of the polynomials Aij = Ylik^ikUjk — ^ij- 
So, U is a critical point of F if and only if dF G spanidAij). Now since A^j — Aji, this is 
the same as asking for a symmetric matrix M such that dF — Y^-^ MijdAij. But: 

J2 MijdAij = MijiUikdUjk + UjkdUik) = 2 Y,iMU)ikdUik 

ij ijk Ik 

Thus the critical point condition reads S — 2MU, so the matrix M — SU*' /2 must be 
symmetric. Now the Hessian of F apphed to a vector X — UY , with Y e 0{N), is: 

Hess{F){X) = -Tr{X' ■ SU' ■ X) = -Tr(F* • U'S ■ Y) 

Thus the Hessian of F is positive definite when U^S is positive definite, which is the 
same as saying that U{U^S)U* = SU* is positive definite, and we are done. □ 

The above result gives rise to the following definition: 

Definition 1.4. A square matrix H e Miv(IR*) is called "almost Hadamard" if U — 
H/ \fN is orthogonal, and the following equivalent conditions are satisfied: 

(1) U is a local maximizer of the 1-norm on 0{N). 

(2) Uij ^ 0, and with Sij = sgn(f/jj), we have SU^ > 0. 

If so is the case, we call H "optimal" if U is a maximizer of the 1-norm on 0{N). 

Let Jn be the flat N x N matrix, having all the entries equal to Also, let us call 

"Hadamard equivalence", or just "equivalence", the equivalence relation on the N x N 
matrices coming from permuting the rows and columns, or multiplying them by —1. 

Proposition 1.5. The class of almost Hadamard matrices has the following properties: 

(1) It contains all the Hadamard matrices. 

(2) It contains the matrix = \/N{2Jn ~ In)- 

(3) It is stable under equivalence, tensor products, and transposition. 
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Proof. All the assertions are clear from definitions: 

(1) This follows either from Proposition 1.2, or from the fact that U = H/\fN is 
orthogonal, and SU^ = HU^ = VnIn is positive. 

(2) First, the matrix U = Kn/\/N is orthogonal, because it is symmetric, and: 

= {2Jn - In? = 44 - ^Jn + In = In 
Also, we have 5* = NJ^ — 21jv, so the matrix SU^ is indeed positive: 

^^7* = (NJn - 2U)(2J;v - In) = (N - 2)Jn + 2(U - Jn) 

(3) For a tensor product of almost Hadamard matrices H = H' ® H" we have U = 
U' ® U" and S = S' ® S", so that U is unitary and SU^ is positive, as claimed. As for the 
assertions regarding equivalence and transposition, these are clear from definitions. □ 

Regarding now the optimal case, we have the following result, from [S]: 

Proposition 1.6. The optimal almost Hadamard matrices are as follows: 

(1) At any N where HC holds, these are the N x N Hadamard matrices. 

(2) At N = 3, these are precisely = \/3{2J^ — I3) and its conjugates. 

Proof. The assertion (1) is clear from Proposition 1.2. For (2) we must prove that for 
U G 0(3) we have \\U\\i < 5, with equality when U is conjugate to U3 = 2J3 — I3. But 
here we can assume U G 50(3), and use the Euler-Rodrigues formula: 

fx^ + y'^ - -t^ 2{yz - xt) 2{xz + yt) \ 

U=\ 2{xt + yz) x^ + z^-y^-t^ 2{zt - xy) 

\ 2{yt-xz) 2{xy + zt) x'^ + t^ - y^ - z'^ J 

Here {x,y,z,t) G comes from the standard cover map ~ SU{2) — > 50(3). Now 
by linearizing, we must prove that for any (x, 2, t) G M'' we have: 

\\U\\,<h{x^ + y^ + z^ + t^) 

The proof of this latter inequality is routine, and the equality situation turns to hold 
indeed exactly for the matrix f/3 = 2J3 — I3 and its conjugates. See [5]. □ 

Finally, let us mention that a version of Proposition 1.2 above, using the Holder in- 
equality, shows that the matrices of type U = H/\fN with H G Mjv(±l) Hadamard 
maximize the p-norm on 0{N) at p G [1, 2), and minimize it at p G (2, 00]. See |5]. Part 
of the above p = 1 results extend to the general setting p G [1, 00] — {2}, and in particular 
to the exponents p = 4 and p = 00, which are of particular interest in connection with 
several quantum physics questions. This will be discussed in a forthcoming paper. 
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2. The circulant case 

In this section we study the almost Hadamard matrices which are circulant. We recall 
that a matrix H e Miq{C) is called circulant if it is of the form: 

7Af-i\ 

lN-2 

70 / 

Let F G U{N) be the Fourier matrix, given by Fij = uj''^ /^/N, where u — e^'^'^l^ . Given 
a vector a G C", we associate to it the diagonal matrix a' = diag{ao, . . . , a^v-i)- 
We will make a heavy use of the following well-known result: 

Proposition 2.1. For a matrix H G Mn{C), the following are equivalent: 

(1) H is circulant, i.e. H^j = jj-i, for a certain vector 7 G C"^. 

(2) H is Fourier-diagonal, i.e. H — FDF* , with D G MjsiifC) diagonal. 

In addition, if so is the case, then with D = \/Na' we have 7 = Fa. 
Proof. (1) =^ (2) The matrix D = F*HF is indeed diagonal, given by: 

kl r 

(2) =^ (1) The matrix H = FDF* is indeed circulant, given by: 

Hij — ^2 FikDkkFjk — — ^y^^u^'^"^''^ Dkk 

k k 

Finally, the last assertion is clear from the above formula of Hij. □ 

Let us investigate now the circulant orthogonal matrices. We let the matrix indices i,j 
vary modulo A^. We denote by T the unit circle in the complex plane. 

Lemma 2.2. For a matrix U G M^ifC), the following are equivalent: 

(1) U is orthogonal and circulant. 

(2) U = Fa'F* with a G satisfying ai = a-i for any i. 

Proof. We will use many times the fact that given a G C^, the vector 7 = Fa is real if 
and only if di = a-i for any i. This foUows indeed from Fa = Fa, with ai = a-i. 

(1) =^ (2) Write Hij = 7j_j with 7 G M^. By using Proposition 2.1 we obtain 
H = FDF* with D = VN a' and 7 = Fa. Now since U = Fa'F* is unitary, so is a', so 
we must have a G T^. Finally, since 7 is real we have ai = a-i, and we are done. 

(2) =^ (1) We know from Proposition 2.1 that U is circulant. Also, from a G we 
obtain that a' is unitary, and so must be U. Finally, since we have ai = a-i, the vector 
7 = Fa is real, and hence we have U G Mjv(M), which finishes the proof. □ 



H 



( 7o 7i 

lN-\ 7o 

V 71 72 
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Let US discuss now the almost Hadamard case. First, in the usual Hadamard case, the 
known examples and the corresponding a-vectors are as follows: 



Proposition 2.3. The known circulant Hadamard matrices, namely 

/-I /ill 



± 



v 



± 



± 



I- 111 11-11 

II- 11 ' 111-1 
1 1 1-1/ \-i iiiy \i-iiiy 

come from the vectors a — ±(1, —1, —1, —1), ±(1, — i, 1, i), ±(1, 1, —1, 1), ±(1, i, 1, —i). 



111-1 

-1111 
1-111 



-1\ 



v 



-1111 

I- 111 

II- 11 



Proof. At N — A the conjugate of the Fourier matrix is given by: 

/II 1 1\ 



F* — 



-1 



□ 



1 -i 
1-11 
\l i -I 

Thus the vectors a — F*j are indeed those in the statement 

We have the following "almost Hadamard" generahzation of the above matrices: 

Proposition 2.4. If — 1 then the vector a — ±(1, — g, — g^, . . . , —q^~^) produces an 
almost Hadamard matrix, which is equivalent to Kn — -\/]V(2Jjv — Iat)- 

Proof. Observe first that these matrices generahze those in Proposition 2.3. Indeed, at 

N = 4 the choices for q are 1, i, —1, —i, and this gives the above a-vectors. 
Assume that the ± sign in the statement is +. With q = a;'", we have: 



N-l 



J2 ^'^<^k 



N-l 



(i+r)k 



N-l 

E 



{i+r)k 



k=0 k=l k=0 

In terms of the standard long cycle {CN)ij — ^i+ij, we obtain: 

H = ViV(2 Jjv - C/) 
Thus H is equivalent to K^, and by Proposition 1.5, it is almost Hadamard. 



□ 



In general, the construction of circulant almost Hadamard matrices is quite a tricky 
problem. At the abstract level, we have the following technical result: 

Lemma 2.5. A circulant matrix H G Mjv(]R*), written H^j = 'jj-i, is almost Hadamard 
if and only if the following conditions are satisfied: 

(1) The vector a — F*7 satisfies a e T-'^. 

(2) With e — sgn(7), pi — ^j.£r^i+r o-nd v — F*p, we have v > 

In addition, if so is the case, then OLi — a^i, pi — p^i and — v_i for any i. 



ALMOST HADAMARD MATRICES 



11 



Proof. According to Definition 1.4 our matrix H is almost Hadamard if any only if the 
matrix U — H/y/N is orthogonal and SU^ > 0, where Sij — sgn{Uij). By Lemma 2.2 
the orthogonality of U is equivalent to the condition (1). Regarding now the condition 
SU^ > 0, this is equivalent to S^U > 0. But, with k = i — r, we have: 



Thus S*U is circulant, with p/y/N as first row. From Proposition 2.1 we get S^U — 
FLF* with L = v' and v = F*p, so S^U > iff z/ > 0, which is the condition (2). 

Finally, the assertions about a, v follow from the fact that Fa, Fv are real. As for the 
assertion about p, this follows from the fact that S^U is symmetric. □ 

Theorem 2.6. For N odd the following matrix is almost Hadamard, 



( 



- COS ^ cos"^ ^ 



cos" 



-1 (Af-l)7r 
N 



— COS 



-1 27r 



-1 Tr_ 
N 



-1 Stt 



COS 
— COS 



1 (iV-l)7r \ 

N 

_1 (Af-2)7r 
N 



\ -COS ' § COS ^ ^ - COS" ^ 



and comes from an a-vector having all entries equal to 1 or —1. 
Proof. Write = 2n + 1, and consider the following vector: 



[-l)n+i for i = 0,l,...,n 

for i = n + l,...,2n 



Let us first prove that {L]^)ij = ''jj-i, where 7 = Fa. With u = e^''*/^ we have: 



2n 



j=0 



j=0 j=l 
N 



Now since N is odd, and since a; = 1, we obtain: 



EM) 

j=0 j=l j=-n 

By computing the sum on the right, with ^ = q^^/^ -^ve get, as claimed: 

2cc;~™ 2^-2"^ 2^-^' 



I'K 



I + UJ' 1 + ^■ 



2i 



- = (-l)^cos-^ 
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In order to prove now that Ljv is almost Hadamard, we use Lemma 2.5. Since the sign 
vector is simply e — (— 1)"'q;, the vector pi — ^^.s^Ti+r is given by: 

2n n n 2n 

,rj 



r=0 3=—n j=—n r=0 

Now since the last sum on the right is {\/NFa)j — Vn^j, we obtain: 

n ^ n n 

^ k=—n 



j=-n 



j=-n 



Thus we have the following formula: 

■1) 



j=—n k=—n 

Let us compute now the vector u = F*p. We have: 

2n 

i 



i=0 j=—n k=—n i=0 



j=—n k=—n 

The sum on the right is NSji, with both j,l taken modulo N, so it is equal to NSjl, 
where L = / for / < n, and L = I — N ior I > n. We get: 



= (-irE(-i) 



L+k, Lk 



A;=— n 



k=—n 



With ^ = e'^*/^, this gives the following formula: 



[-ly 



-11 



2u 



-nL 



In terms of the variable ^ = e^^^^ , we obtain: 



1 + 



(-i)V 



_1 -^TT 



r /-r COS 



Now since L e [— n, n], all the entries of u are positive, and we are done. 



□ 



At the level of examples now, at A?" = 3 we obtain the matrix L3 = —K^. At N — 5 we 
obtain a matrix having as entries 1 and x — — cos~^ t, y — cos~^ ^: 



(1 X y y x\ 

X 1 X y y 

y X 1 X y 

y y X 1 X 

\x y y X Ij 
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Let US look now more in detail at the vectors a G T appearing in Proposition 2.4 and 
in the proof of Theorem 2.6. In both cases we have = u^^ for a certain r G N. We 
conjecture that all the circulant almost Hadamard matrices appear in this way: 

Conjecture 2.7. Any circulant almost Hadamard matrix must come from a vector a G 
having the property that o? is formed by roots of unity in a progression. Moreover, 
the passage o? ^ a is via the square root operation with respect to a half-plane. 

It would be nice to further refine this statement, as to have a complete list of possible 
circulant almost Hadamard matrices. But we have no further results here. 



In this section we study the almost Hadamard matrices having only two entries, H G 
M]\f{x,y), with x,y eR. As a first remark, the usual Hadamard matrices H G M]y{±l) 
are of this form. However, when trying to build a combinatorial hierarchy of the two-entry 
almost Hadamard matrices, the usual Hadamard matrices stand on top, and there is not 
so much general theory that can be developed, as to cover them. 

We will therefore restrict attention to the following special type of matrices: 

Definition 3.1. An {a,b,c) pattern is a matrix M G Mp^^x.y), with N = a + 2b + c, such 
that, in any two rows, the number of x/y/x/y sitting below x/x/y/y is a/b/b/c. 

In other words, given any two rows of our matrix, we are asking for the existence of a 
permutation of the columns such that these two rows become: 

X . . .X X . . .X y . . .y y ■ ■ - y 
x^^^^ y...y y---y 

a I, b c 

Oberve that the Hadamard matrices do not come in general from patterns. However, 
there are many interesting examples of patterns coming from block designs [H], |35] : 

Definition 3.2. A {v,k,X) symmetric balanced incomplete block design is a collection B 
of subsets of a set X , called blocks, with the following properties: 

(1) |X| = \B\ = V. 

(2) Each block contains exactly k points from X . 

(3) Each pair of distinct points is contained in exactly A blocks of B. 

The incidence matrix of a such block design is the f x t> matrix defined by: 



3. The two-entry case 




The connection between designs and patterns comes from: 
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Proposition 3.3. // = a + 2b + c then the adjacency matrix of any {N, a + b,a) 
symmetric balanced incomplete block design is an {a, b, c) pattern. 

Proof. Indeed, let us replace the — 1 values in the adjacency matrix M by abstract x — y 
values. Then each row of M contains a + b copies of x and b + c copies of y, and since 
every pair of distinct blocks intersect in exactly a points, cf. we see that every pair of 
rows has exactly a variables x in matching positions, so that M is an (a, b, c) pattern. □ 




Figure 1. The Fano plane 

As a first example, consider the Fano plane. The sets X, B of points and lines form a 
(7, 3, 1) block design, corresponding to the following (1, 2, 2) pattern: 

y\ 



(x 

X 
X 

y 
y 
y 

\y 



X 

y 
y 

X 
X 

y 
y 



X 

y 
y 
y 
y 

X 
X 



y 

X 

y 

X 

y 

X 

y 



y 

X 

y 
y 

X 

y 

X 



y 
y 

X 
X 

y 
y 

X 



y 

X 

y 

X 
X 

y) 



Now remember that the Fano plane is the projective plane over F2 = {0, 1}. The same 
method works with F2 replaced by an arbitrary finite field F^, and we get: 

Proposition 3.4. Assume that q = p^ is a prime power. Then the point-line incidence 
matrix of the projective plane over Fg is a (1, g, — q) pattern. 

Proof. The sets X, B of points and lines of the projective plane over F^ are known to form 
a (g^ + g + 1, g + 1, 1) block design, and this gives the result. □ 

There are many other interesting examples of symmetric balanced incomplete block 
designs, all giving rise to patterns, via Proposition 3.3. For instance the famous Paley 
biplane [11], pictured below, is a (11, 5, 2) block design, and hence gives rise to a (2, 3, 3) 
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pattern. When assigning certain special values to the parameters y we obtain a 11 x 11 
almost Hadamard matrix, that we believe to be optimal. See section 4 below. 




Figure 2. The Paley biplane 

We consider now the problem of associating real values to the symbols x, y in an (a, b, c) 
pattern such that the resulting matrix U{x,y) is orthogonal. 

Lemma 3.5. Given a,b,c & N, there exists an orthogonal matrix having pattern (a, b, c) 
iffb"^ > oc. In this case the solutions are U{x,y) and —U{x,y), where 

_ t _ 1 

""''Tkt + Yy ^~VHt + i) 

where t = (b± \Jb'^ — ac) / a can be any of the solutions of at^ — 2bt + c = 0. 

Proof. First, in order for f/ to be orthogonal, the following conditions must be satisfied: 

ax"^ + 2bxy + cy'^ = 0, (a + b)x'^ + {b + c)y'^ = 1 

The first condition, coming from the orthogonality of rows, tells us that t = —x/y must 
be the variable in the statement. As for the second condition, this becomes: 

2_ 1 _ 1 

^ ~ (a + 6)^2 + (6 + c) ~ (at2 + c) + (6^2 + h) ~ 2bt + bt^ + b ~ HtTTf 

This gives the above formula of y, and hence the formula of a; = —ty as well. □ 

Lemma 3.6. Let U = U{x,y) be orthogonal, corresponding to an (a, 6, c) pattern. Then 
H = ^/NU is almost Hadamard iff {N{a - b) + 2b)\x\ + {N{c - b) + 2b)\y\ > 0. 

Proof. We use the criterion in Definition 1.4 (2). So, let Sij = sgn{Uij). Since any row of 
U consists of a + 6 copies of x and b + c copies of y, we have: 

{SU% = J2 sgn{Uik)Uik = (a + b)\x\ + {b + c)\y\ 

k 
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Regarding now {SU^)ij with i ^ j, we can assume in the computation that the i-th and 
j-th row of U are exactly those pictured after Definition 3.1 above. Thus: 

{SU% = J]sgn([/,fc)^7,fc 

k 

— a sgn(a;)a; + b sgn{x)y + b sgn{y)x + c sgn(y)y 
= a\x\ — b\y\ — b\x\ + c\y\ 
= {a-b)\x\ + {c-b)\y\ 

We obtain the following formula for the matrix SU* itself: 

SU' = 2b{\x\ + \y\)lM + {{a-b)\x\ + {c-b)\y\)NJN 

= 2b{\x\ + |y|)(ljv - Jn) + {{N{a - b) + 2b)\x\ + {N{c - b) + 26)|y|)) Jjv 

Now since the matrices In — Jn, Jn are orthogonal projections, we have SU* > if and 
only if the coefficients of these matrices in the above expression are both positive. Since 
the coefficient of Ijv — ^Af is clearly positive, the condition left is: 

{N{a -b)+ 2b)\x\ + {N{c - b) + 2b)\y\ > 

So, we have obtained the condition in the statement, and we are done. □ 

Proposition 3.7. Assume that a,b,c satisfy c> a and b{b — 1) = ac, and consider 
the {a,b,c) pattern U — U{x,y), where: 



a 



+ {l-a-b)\/b b+{a + b)Vb 



Na ' ^ Nb 

Then H — y/NU is an almost Hadamard matrix. 

Proof. We have b'^ — ac — b, so Lemma 2.5 applies, and shows that with t — (b — Vb)/a 
we have an orthogonal matrix U — U{x,y), where: 

t 1 
^ = 777. 77' 



Vb{t + iy Vb{t + 1) 

In order to compute these variables, we use the following formula: 

(a + b)'^ - b^a^ + b'^ + 2ab-b^a^ + 2ab + ac^ Na 
This gives indeed the formula of y in the statement: 

a {a + b)\/b + b 

^~ {a + b)Vb-b ~ Nb 
As for the formula of x, we can obtain it as follows: 

_ (V^ - b){{a + b)Vb + b) _ a + {l-a- b)Vb 
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Let us compute now the quantity appearing in Lemma 3.6. We have: 

N{a-b) + 2b = {a + 2b + c){a - b) + 2b 

= a? + ab- 26^ + ac-bc + 2b 
= + ab — ac — be 
= (a — c)(a + 6) 

Similarly, N{c — b) + 2b = (c — a)(c + 6), so the quantity in Lemma 3.6 is Ky, with: 

K = {a-c){a + b)t + {c-a){c + b) 
= (c- a){c + b- {a + b)t) 

{ac + ab - {a + b){b - Vb)) 
{{ac-b^) + {a + b)Vb) 
iia + b)Vb-b) 

Since this quantity is positive, Lemma 3.6 applies and gives the result. □ 

Theorem 3.8. Assume that q = is a prime power. Then the matrix In G MN{x,y), 
where N = + q + 1 and 

l-?v^ ? + + 

^ = — 7^^' y ~ 



c 




a 




a 




c 




a 




a 




c 




a 


a 



'N qVN 
having (l,g,g^ — q) pattern coming from the point-line incidence of the projective plane 



over ¥q is an almost Hadamard matrix. 

Proof. Indeed, the conditions c > a and b{b — 1) = ac needed in Proposition 3.7 are 
satisfied, and the variables constructed there are x' = x/\/N and y' = y/\fN. □ 

There are of course many other interesting examples of two-entry almost Hadamard 
matrices, all worth investigating in detail, but we will stop here. Indeed, the main purpose 
of the reminder of this paper is to provide a list of almost Hadamard matrices which are 
"as optimal as possible", at = 2, 3, . . . , 13, and in order to establish this list, we will 
just need the incidence matrices Jat, plus the matrix coming from the Paley biplane. 

Let us mention however two more important aspects of the general theory: 

(1) The series 1^ is the particular case of a 2-parameter series Indeed, associated 
to q = p^ and G N is a certain {[d + 2]g, [d + 1]^, [d]g) block design coming from 
Fg, where [e]g = {q^ — l)/(g — 1), cf. [H], [35]. Thus by Proposition 3.7 we obtain 
an almost Hadamard matrix having pattern (^fj-, 5''^, q'^{q — !))• 

(2) Trying to find block designs and patterns leads to the following chess problem: 
consider a. N x N chessboard, take NM rooks with M < N/2, fix an integer 
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K < M, and try to place all the rooks on the board such that: (a) there are 
exactly M rooks on each row and each column of the board, and (b) for any pair 
of rows or columns, there are exactly K pairs of mutually attacking rooks. 

Finally, let us mention that there are many questions raised by the almost Hadamard 
matrices, at the quantum algebraic level. For instance the symmetries of Hadamard matri- 
ces are known to be described by quantum permutations [3j, and it would be interesting to 
have a similar result for the almost Hadamard matrices. This might probably bring some 
new ideas on the "homogeneous implies quantum homogeneous" question, raised in [2], 
and having connections with the finite projective planes. Also, an interesting link between 
mutually unbiased bases, complex Hadamard matrices and affine planes was emphasized 
in [H], but its relation with our present investigations is not known yet. 

4. List of examples 

In this section we present a list of examples of almost Hadamard matrices, for small 
values of N. Since we are mainly interested in the optimal case, our examples will be 
chosen to be "as optimal as possible", i.e. will be chosen as to have big 1-norms. 

So, let us first compute the 1-norms for the examples that we have. In what follows 
Ltv is the matrix found in Theorem 2.6, and 1^ is the matrix found in Theorem 3.8. 

Theorem 4.1. The 1-norms of the basic examples of almost Hadamard matrices are: 

(1) Hadamard case: | |if/-\/iV| |i = N\/l^, for any H G Mn{±1) Hadamard. 

(2) Basic series case: \\Kn/\/N\\i = 3N — A, where = ^/N{2JN — In)- 

(3) Circulant series case: \\Ln/\/N\\i = ^NlogN + 0{N). 

(4) Incidence series case: | |/Ar/-\/iV| |i = (g^ — g — 1) + 2g(g + l)y/g. 

Proof. The first two assertions are clear. For the third one, with = 2r;, + 1 we have: 
\\L^/VN\h = 2 X: cos-^ ^ + 0(iV) = 2 X: sin-^ ^^^^ + 0{N) 

1=0 i=0 

Now by using sinx ~ x and Yli=i 1/"^ = log A; + 0(1) we obtain, as claimed: 
\\Ln /v^Wi = — E TTTT + 0(iV) = — log AT + 0{N) 

1=0 

As for the last assertion, let first U be the matrix in Lemma 3.6. We have: 

\\U\U = N{{a + b)\x\ + {b + c)\y\) 
In the particular case of the orthogonal matrices in Proposition 3.7, we get: 

N {a + b){2a + 2b-2) r- 
\\U\\i = (c-a) + ^ '-Vb 

Now with a = 1, b = q, c = — q, this gives the formula in the statement. □ 
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Observe that at big the above matrices Kn, Ln, In are far from being optimaL In 
fact, with N ^ oo, the corresponding orthogonal matrices don't even match the average 
of the 1-norm on 0{N), which, according to |5], is ~ cN\/N, with c = 0.797.. 

With these ingredients in hand, let us discuss now the various examples: 

At = 2 we have the Walsh matrix H2, which is of course optimal. 

At = 3 we have the almost Hadamard matrix K^, which is optimal: 




At N = 4 we have the Hadamard matrix H4 ~ K4, once again optimal. 
At = 5 we have the matrix K^, that we believe to be optimal as well: 



/-3 2 2 2 2 \ 

2 -3 2 2 2 

2 2-3 2 2 

2 2 2 -3 2 

2 2 2 2 -3/ 



At A^ = 6 it is plausible that the optimal AHM is simply H2: 



f-l 2 2 -1 2 2 \ 

2-12 2-12 

2 2 -1 2 2 -1 

-12 2 1-2-2 
2-12-21-2 
\ 2 2-1-2-21/ 



At A^ = 7 we have the incidence matrix of the Fano plane {x = 2 — 4-\/2, y = 2 + 3\/2): 



2V7 



fx X X y y y y\ 

X y y X X y y 

X y y y y X X 

y X y X y y X 

y X y y X X y 

y y X y X y X 

\y y X X y X yj 



At A^ = 8 we have the third Walsh matrix = H2® i^4, of course optimal. 

At A^ = 9 we just have the matrix ® K3, which can be shown not to be optimal. 
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At = 10 we believe that the matrix H2 is optimal: 



AT. ®Ho = ^ 



1 



/-3 
2 

2 

2 

2 

-3 
2 
2 
2 
V2 



2 
2 
2 
2 

-3 
2 
2 
2 



2 
2 

-3 
2 
2 
2 
2 

-3 
2 
2 



2 
2 
2 

-3 
2 
2 
2 
2 



2 
2 
2 
2 

-3 
2 
2 
2 
2 



-3 2 



2 

2 
2 
2 
3 



2 
2 
2 



-2 
-2 
-2 

-2 



-2 
-2 

-2 



2 
2 

-3 
2 
2 

-2 
-2 

3 

-2 

-2 



2 
2 
2 



-2 
-2 
-2 



2 \ 
2 

2 

2 

-3 
-2 
-2 
-2 
-2 
3/ 



At A" = 11 we have the matrix of the Paley biplane {x — & — 12\/3, y — & + 10-\/3): 



11 
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y 


y 


X 


y 


y 


X 


X 


X 


y 


y 


y 


X 


y 


y 


X 


\x 


y 


X 


X 


X 


y 


y 


y 


X 


y 


y > 



At A^ = 12 wc have the Sylvester Hadamard matrix 5*12, of 
At A^ = 13 we have the incidence matrix of P(F3) (x = 3 — 



course optimal. 
9^3, y = 3 + 4^3): 



'13 



3V13 
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y 


X 


y 


y J 
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N 


matrix 


1-norm (formula) 


1-norm (numeric) 


— 

Ny/N 


remarks 


2 


H2 


2V2 


2.828 


2.828 


Hadamard 


3 




5 


5.000 


5.196 


optimal 


4 


Ki 


8 


8.000 


8.000 


Hadamard 


5 


K, 


11 


11.000 


11.118 




6 


K30H2 


10v^8 


14.142 


14.697 




7 


h 


1 + I2V2 


17.971 


18.520 




8 


Hs 


16^2 


22.627 


22.627 


Hadamard 


9 


? 




> 26.513 


27.000 




10 




22 


31.113 


31.623 




11 


Pu 


1 + 20^3 


35.641 


36.483 




12 




24^3 


41.569 


41.569 


Hadamard 


13 


/l3 


5 + 24^3 


46.569 


46.872 





Table 1. Almost Hadamard matrices H G Mn{W), chosen as for the corre- 
sponding orthogonal matrices U = H/ \/N to have big 1-norm. All matrices 
are believed to be optimal. The lower bound for the maximum of the 1-norm 
on 0(9), which is not optimal, was obtained by numerical simulation. 



The norms of the above matrices can be computed by using the various formulae in 
Theorem 4.1 and its proof, and the results are summarized in Table 1. 

5. Conclusion 

We have seen in this paper that the Hadamard matrices are quite nicely generalized by 
the almost Hadamard matrices (AHM), which exist at any given order G N. Our study 
of these matrices, which was for the most of algebraic nature, turns to be related to several 
interesting combinatorial problems, notably to the Circulant Hadamard Conjecture. 

We believe that the AHM can be used as well in connection with several problems in 
quantum physics, in a way somehow similar to the way the complex Hadamard matrices 
(CHM) are used. Indeed, since the CHM exist as well at any given order N E N, these 
matrices proved to be useful in several branches of quantum physics. For instance in 
quantum optics they are sometimes called the "Zeilinger matrices" , as they can be applied 
to design symmetric linear multiports, used to split the beam into A^ parts of the same 
intensity and to analyze interference effect [33] , [22] . In the theory of quantum information 
one uses quantum Hadamard matrices to construct mutually unbiased bases (MUB) [21] . 
[16] to design teleportation and dense coding schemes. As shown in the seminal work of 
Werner [38] these two problems are in fact equivalent and also equivalent to construction 
of unitary depolarisers and maximally entangled bases [39] . 
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Although quantum mechanics in a natural way relays on a complex Hilbert space it 
is often convenient to study a simplified problem and restrict attention to the subset of 
real quantum states only. Such an approach can be useful in theoretical investigations 
of quantum entanglement ^13] or also in experimental studies on engineering of quantum 
states, as creating of a real state by an orthogonal rotation usually requires less effort 
than construction of an arbitrary complex state. 

There exists therefore a natural motivation to ask similar problems concerning e.g. 
unbiased bases and teleportation schemes in the real setup. For instance, it is known that 
for any N there exist < N/2 + 1 real MUB and for most dimensions their actual number 
is not larger than 3, cf. [IQl, while for any prime N there exist + 1 complex MUB. Note 
that the real MUB and Hadamard matrices are closely related to several combinatorial 
problems [19], [29], [30]. In the case the maximal number of real MUB does not exist one 
can search for an optimal set of real bases which are approximately unbiased. 

Let us now return to the construction of symmetric multiports which relay on the com- 
plex Hadamard matrices [33] , [22] . In the case of a fc-qubit system there exist Hadamard 
matrices of order = 2^, so one may use real orthogonal matrices for this purpose. 
However, already for = 3, 5, 6, 7 real Hadamard matrices do not exist, so in these di- 
mensions, if one restricts the rotations to orthogonal matrices, there are no real symmetric 
multiports. Therefore, for these dimensions one may always raise the following question: 
what is the optimal, approximate solution of the problem, if one is allowed to use only 
real states and orthogonal matrices? The almost Hadamard matrices analyzed in this 
paper are directly applicable for such a class of problems. 

References 

[1] A. Alvarez, J. A. Armario, M.D. Frau, E. Martin and A. Osuna, Error correcting codes from quasi- 
Hadamard matrices, Lecture Notes in Comput. Sci. 4547 (2007), 294-302. 

[2] T. Banica, Quantum automorphism groups of homogeneous graphs, J. Funct. Anal. 224 (2005), 
243-280. 

[3] T. Banica, Quantum permutations, Hadamard matrices, and the search for matrix models, arxiv: 
1109.4888. 

[4] T. Banica, J. Bichon and J.-M. Schlenker, Representations of quantum permutation algebras, /. 
Funct. Anal. 257 (2009), 2864-2910. 

[5] T. Banica, B. Collins and J.-M. Schlenker, On orthogonal matrices maximizing the 1-norm, Indiana 
Univ. Math. J. 59 (2010), 839-856. 

[6] T. Banica, B. Collins and J.-M. Schlenker, On polynomial integrals over the orthogonal group, J. 
Comhin. Theory Ser. A 118 (2011), 778-795. 

[7] T. Banica and J.-M. Schlenker, Combinatorial aspects of orthogonal group integrals, Internat. J. 
Math. 22 (2011), 1611-1646. 

[8] I. Bengtsson, W. Bruzda, A. Ericsson, J. -A. Larsson, W. Tadej and K. Zyczkowski, Mutually unbi- 
ased bases and Hadamard matrices of order six, J. Math. Phys. 48 (2007), 1-33. 

[9] I. Bengtsson and A. Ericsson, Mutually unbiased bases and the complementarity polytope. Open 
Syst. Inf. Dyn. 12 (2005), 107-120. 



ALMOST HADAMARD MATRICES 



23 



[10] P.O. Boykin, M. Sitharam, M. Tarifi and P. Wocjan, Real mutually unbiased bases, arxiv: 
quant -ph/0502024. 

[11] E. Brown, The fabulous (11,5,2) biplane. Math. Mag. 77 (2004), 87-100. 

[12] A.T. Butson, Generalized Hadamard matrices, Proc. Amer. Math. Soc. 13 (1962), 894-898. 

[13] CM. Caves, C.A. Fuchs and P. Rungta, Entanglement of formation of an arbitrary state of two 
rebits. Found. Phys. Lett. 14 (2001), 199. 

[14] C.J. Colbourn and J.H. Dinitz, Handbook of combinatorial designs, CRC Press (2007). 

[15] B. Collins and P. Sniady, Integration with respect to the Haar measure on unitary, orthogonal and 
symplectic groups. Comm. Math. Phys. 264 (2006), 773-795. 

[16] T. Durt, B.-G. Englert, I. Bengtsson and K. Zyczkowski, On mutually unbiased bases, Int. J. Quan- 
tum Information 8 (2010), 535-640. 

[17] U. Haagerup, Orthogonal maximal abelian *-subalgebras of the nxn matrices, in "Operator algebras 
and quantum field theory". International Press (1997), 296-323. 

[18] J. Hadamard, Resolution d'une question relative aux determinants. Bull. Set. Math. 2 (1893), 240- 
246. 

[19] W.H. Holzman, H. Kharaghani and W. Orrick, On real unbiased Hadamard matrices, Contemp. 

Math. 531 (2010), 243 250. 
[20] K.J. Horadam, Hadamard matrices and their applications, Princeton Univ. Press (2007). 
[21] I.D. Ivanovic, Geometrical description of quantal state determination, J. Phys. A 14 (1981), 3241- 

3245. 

[22] I. Jex, S. Stenholm and A. Zeilinger, Hamiltonian theory of a symmetric multiport. Opt. Com,mun. 
117 (1995), 95-101. 

[23] V.F.R. Jones, On knot invariants related to some statistical mechanical models. Pacific J. Math. 
137 (1989), 311-334. 

[24] H. Kharaghani and B. Tayfeh-Rezaie, A Hadamard matrix of order 428, J. Comhin. Des. 13 (2005), 
435-440. 

[25] J. Konvalina and R.H. Kosloski, Cyclic quasi- Hadamard matrices, Utilitas Math. 12 (1977), 273-287. 
[26] T.Y. Lam and K.H. Leung, On vanishing sums of roots of unity, J. Algebra 224 (2000), 91-109. 
[27] W. de Launey, On the non-existence of generalized weighing matrices, Ars Comhin. 17 (1984), 
117-132. 

[28] W. de Launey and D.M. Gordon, A comment on the Hadamard conjecture, J. Comhin. Theory Ser. 
A 95 (2001), 180 184. 

[29] N. LeCompte, W.J. Martin and W. Owens, On the equivalence between real mutually unbiased bases 
and certain class of association schemes, European J. Comhin. 31 (2010), 1499-1512. 

[30] M. Matolcsi, I.Z. Ruzsa and M. Weiner, Real and complex unbiased Hadamard matrices, arxiv: 
1201.0631. 

[31] K.-H. Park and H.-Y. Song, Quasi-Hadamard matrices, Proc. ISIT 2010, Austin, TX (2010). 
[32] S. Popa, Orthogonal pairs of *-subalgebras in finite von Neumann algebras, J. Operator Theory 9 
(1983), 253-268. 

[33] M. Reck, A. Zeilinger, H.J. Bernstein and P. Bertani, Experimental realization of any discrete unitary 

operator, Phys. Rev. Lett. 73 (1994), 58-61. 
[34] H.J. Ryser, Combinatorial mathematics, Wiley (1963). 

[35] D.R. Stinson, Combinatorial designs: constructions and analysis. Springer- Vcrlag (2006). 

[36] J.J. Sylvester, Thoughts on inverse orthogonal matrices, simultaneous sign-successions, and tesse- 

lated pavements in two or more colours, with applications to Newton's rule, ornamental tile-work, 

and the theory of numbers, Phil. Mag. 34 (1867), 461-475. 



24 



TEODOR BANICA, ION NECHITA, AND KAROL ZYCZKOWSKI 



[37] W. Tadej and K. Zyczkowski, A concise guide to complex Hadamard matrices, Open Syst. Inf. Dyn. 
13 (2006), 133-177. 

[38] R.F. Werner, All teleportation and dense coding schemes, J. Phys. A 34 (2001), 7081-7094. 
[39] A. Wojcik, A. Grudka and R.W. Chhajlany, Generation of inequivalent generalized Bell bases. Quan- 
tum Inf. Process. 2 (2003), 201-206. 

T.B.: Department of Mathematics, Cergy-Pontoise University, 95000 Cergy-Pontoise, 
France, teodor.bcinica@u-cergy.fr 

I.N.: CNRS, Laboratoire de Physique Theorique, IRSAMC, Universite de Toulouse, 
UPS, 31062 Toulouse, France. nechitaOirsamc.ups-tlse.fr 

K.Z.: Institute of Physics, Jagiellonian University, Cracow and Center for Theoret- 
ical Physics, Polish Academy of Sciences, Warsaw, Poland. karolQtatry. if .uj .edu.pl 



